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Photoperiod sensitivity is an important trait related to crop adaptation and ecological breeding in common 
buckwheat (Fagopyrum esculentum Moench). Although photoperiod sensitivity in this species is thought to 
be controlled by quantitative trait loci (QTLs), no genes or regions related to photoperiod sensitivity had been 
identified until now. Here, we identified QTLs controlling photoperiod sensitivity by QTL analysis in a seg- 
regating F 4 population («=100) derived from a cross of two autogamous lines, 02AL1 13(Kyukei 
SC2)LH.self and C0408-0 RP. The F 4 progenies were genotyped with three markers for photoperiod- 
sensitivity candidate genes, which were identified based on homology to photoperiod-sensitivity genes in 
Arabidopsis and 76 expressed sequence tag markers. Among the three photoperiod-sensitivity candidate genes 
(FeCCAl, FeELFS and FeCOL3) identified in common buckwheat, FeELF3 was associated with photoperiod 
sensitivity. Two EST regions, Fest_L0606_4 and Fest_L0337_6, were associated with photoperiod sensitivity 
and explained 20.0% and 14.2% of the phenotypic variation, respectively. For both EST regions, the allele 
from 02AL113(Kyukei SC2)LH. self led to early flowering. An epistatic interaction was also confirmed be- 
tween Fest_L0606_4 and Fest_L0337_6. These results demonstrate that photoperiod sensitivity in common 
buckwheat is controlled by a pathway consisting of photoperiod-sensitivity candidate genes as well as mul- 
tiple gene action. 

Key Words: common buckwheat, EST markers, photoperiod-sensitivity candidate genes, photoperiod sensi- 
tivity, QTL analysis. 



Introduction 

Common buckwheat (Fagopyrum esculentum Moench, 
2n = 16), a pseudocereal, has been widely cultivated in Asia, 
Europe and North America. In Japan, common buckwheat 
has been cultivated throughout most of the country, with cul- 
tivars classified into three ecotypes: summer, intermediate, 
and late-summer. The late-summer cultivars are cultivated 
in low-latitude regions where photoperiod is around 12 to 
14 h during the cropping season, and their flowering is gen- 
erally sensitive to photoperiod. The summer cultivars are 
cultivated in high-latitude regions where photoperiod is 
around 14.5 to 15.5 h during the cropping season (summer 
crop season), and they are generally insensitive to photo- 
period. Intermediate cultivars show moderate photoperiod 
sensitivity (Namai 1990). Common buckwheat is consid- 
ered to be a short-day plant, with the threshold between flow- 
ering and non-flowering at a photoperiod of about 14.5 h 
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(Nagatomo 1961). Under long-day conditions, (i.e., photo- 
period longer than 14.5 h), late-summer cultivars showed 
wide variation in flowering time, which resulted in low yield 
because of the increased incidence of malformed flowers 
and the decrease in ripening rate in the late-flowering indi- 
viduals (Nagatomo 1961, Nakamura and Nakayama 1950, 
Sugawara 1958). The intermediate and summer cultivars are 
considered to have been derived from the late-summer culti- 
vars through both natural and artificial selection (Matano 
and Ujihara 1979). Minami and Namai (1986a, 1986b) sug- 
gested that summer cultivars differentiated from late- 
summer cultivars through the selection of early-flowering 
plants under long-day conditions; this selection was part of 
the domestication of buckwheat for climatic adaptation to 
the northern part of Japan. 

Iwata et al. (2005) reported that allelic richness (i.e., the 
degree of polymorphism of genetic markers) has decreased 
in Japanese common buckwheat varieties cultivated in high- 
latitude regions during selection for early flowering time, 
and they suggested that the decrease of allelic richness at 
neutral loci was caused by bottleneck effects caused by di- 
rectional selection based on photoperiod sensitivity. Ohsawa 
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(1997) reported that the heritability of flowering period was 
high (95.1%), suggesting that photoperiod sensitivity was a 
key feature of adaptation and differentiation in common 
buckwheat. Photoperiod-sensitivity genes, however, have 
not been identified in common buckwheat. 

Minami (1985) suggested that photoperiod sensitivity in 
common buckwheat was controlled by a number of genes. In 
several crops, photoperiod sensitivity is considered to be in- 
herited quantitatively (Putterill et al. 2004, Yano et al. 
2001). Genetic regulation of flowering induced by photo- 
period has been studied inArabidopsis and rice {Oryza sativa 
L.), and many genes related to photoperiod sensitivity have 
been detected as quantitative trait loci (QTLs) and isolated 
by positional cloning (Corbesier and Coupland 2005, 
Hayama and Coupland 2003). These genes compose the 
photoperiod-sensitivity pathway (Putterill et al. 2004), and 
have been confirmed among many species as orthologous 
genes (Hayama et al. 2002, 2003, Hayama and Coupland 
2003, Kojima et al. 2002, Liu et al. 2001a, Nemoto et al. 
2003, Yano et al. 2000). Moreover, photoperiod sensitivity 
in both long-day and short-day plants appears to be con- 
trolled by a common mechanism (Putterill et al. 2004). 
Based on this information, we hypothesized that photo- 
period sensitivity of common buckwheat would also be con- 
trolled by the action of multiple genes orthologous to those 
of the photoperiod-sensitivity pathways in other species. 

To avoid the complex genetic patterns caused by the allog- 
amous (cross-pollinating) reproductive system of most com- 
mon buckwheat, we used autogamous (self-pollinating) lines 
showing differences in photoperiod sensitivity for simplicity 
of analysis. We used expressed sequence tag (EST) markers 
derived from cDNAs of candidate photoperiod-sensitivity 
genes to detect photosensitivity QTLs most efficiently. 

In the present study, QTLs related to photoperiod sensi- 
tivity in common buckwheat were identified in the following 
three steps: (1) evaluation of variation in photoperiod sensi- 
tivity between two autogamous parental lines, (2) linkage 
map development using EST markers, and (3) QTL analysis 
for photoperiod sensitivity. 

Materials and Methods 

Photoperiod sensitivity of parental lines 

For detecting QTLs related to photoperiod sensitivity, we 
employed two parental lines, 02AL1 13(Kyukei SC2)LH.self 
(KYU) and C0408-0RP (CAN). KYU is an autogamous line 
bred at the National Agricultural Research Center for 
Kyushu Okinawa Region, and CAN is an autogamous line 
bred at Kade Research Ltd. (Canada). In general, differences 
in photoperiod sensitivity between these two lines was visi- 
ble under photoperiods longer than 14.5 h. We used two 
photoperiod conditions, a natural long-day condition (day 
length of -14.5 h) and an artificial long-day condition 
(15.5 h), to detect the difference in photoperiod sensitivity 
between KYU and CAN. For the natural long-day condition, 
60 seeds each of KYU and CAN were sowed by 12 seeds per 



one planter (19x59x 16 cm) and cultivated in growth cabi- 
net without supplemental light at the University of Tsukuba. 
For the artificial long-day condition, 12 seeds each of KYU 
and CAN were sowed in a planter and cultivated in a growth 
cabinet with day length of 15.5 h. Metal halide lamps 
(MLBOC400C-U) (Mitsubishi Electric Corp., Tokyo, 
Japan) were used to control the day length at artificial long- 
day condition. 

The investigation of photoperiod sensitivity was per- 
formed as follows. The dates of cotyledon development and 
of the first flower flowering were recorded for each individ- 
ual. The number of days from the expansion of cotyledons to 
the first flower flowering was defined as days-to-flowering, 
which was used as a measure of the photoperiod sensitivity 
of each individual. The measurements were performed every 
day until 100 days after sowing. At 100 days after sowing, 
individuals that had not yet flowered were classified as 'non- 
flowering' and these were excluded from the following anal- 
yses. To evaluate differences in photoperiod sensitivity be- 
tween the populations, we performed analysis of variance 
(ANOVA) and t tests based on the days-to-flowering under 
each photoperiod condition by using the program IMP 6.0 
(SAS Institute Inc., NC, USA). 

cDNA library construction and search for photoperiod- 
sensitivity candidate gene regions 

We developed sets of 362 cDNA clones derived from in- 
florescences and 1920 cDNA clones derived from leaves. 
362 cDNA clones derived from inflorescences was provided 
by the Niigata University of Pharmacy and Applied Life Sci- 
ences. To develop 1920 cDNA clones derived from leaves, 
we used 'Miyazakizairai' (late-summer cultivar) for plant 
material and sampled leaves six times at 4-h intervals from 
individuals that were either early flowering or non-flowering 
under the 15.5 h day-length condition because it is clarified 
that amount of photoperiod-sensitivity genes expression in 
Arabidopsis was changed over time. Total RNA was extract- 
ed by using RNeasy Plant Mini Kit (Qiagen Inc., CA, USA) 
following manufacturer's instructions. cDNA library was di- 
rectionally constructed (5' SfilA, 3' SfilB) in pDNR-LIB us- 
ing the Creator SMART cDNA Library Construction Kit 
(Clontech Laboratories, Inc., CA, USA) and size-fractionated 
by using CHROMA SPLN+TE-1000 (Clontech Laborato- 
ries, Inc.). Libraries were manually arrayed in 96-well mi- 
crotiter plates. Glycerol stocks of overnight cultures were 
prepared in 96-well format. Plasmid DNAs were extracted 
and BigDye Terminator (Life Technologies Corp., NY, 
USA) cycle sequenced on ABI PRISM 3730x1 Genetic Ana- 
lyzer (Life Technologies Corp.) by using conventional pro- 
cedures and the following primers: M13 forward primer (5'- 
GTAAAACGACGGCCAGT-3') and M13 reverse primer 
(5 '-AAACAGCTATGACC ATGTTC A). Vector sequences 
and low quality regions were trimmed from EST sequences 
manually before removing overlapping regions by using 
BLASTN, and then 863 clones (139 inflorescence-derived 
clones, designated as the Fest_F group, and 724 leaf-derived 
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clones, designated as the FestL group) were used for fur- 
ther analysis. To search for photoperiod-sensitivity candi- 
date gene regions, a homology search against previously re- 
ported sequences to Arabidopsis was performed by using 
tBLASTX. Matches were considered to be significant when 
the smallest sum probability (p) was less than 0.0001 and bit 
scores were greater than 100. 

Development of the mapping population 

We developed a mapping population consisting of 1 00 F 4 
lines derived by single-seed descent from a cross between 
KYU as the female parent and CAN as the male parent. Be- 
cause we observed stunted growth in the F 5 generation, we 
used the F 4 generation as the population for linkage mapping 
and QTL analysis. 

Total DNA samples of both parental lines and the F 4 lines 
were extracted by the following method. About 100 mg leaf 
tissue was mixed with 500 ul lysis buffer (0.3% sodium 
dodecyl sulfate, 20 mM Tris-HCl (pH8.0 at 25°C), 5mM 
EDTA, 400 mM NaCl), and incubated at 65 °C for lOmin 
with 5 ul of lOmg/ml Proteinase K (Wako Pure Chemical 
Ind., Ltd., Osaka, Japan) and lOOmg/ml RNase A (Qiagen 
Inc.). The DNA was purified with chloroform-octanol 
(24 : 1) and precipitated with isopropanol. The pellet was 
washed with 70% ethanol and resuspended in 100 Tris- 
EDTA buffer. 

Marker development and linkage map construction 

Primer sets corresponding to three candidate genes (iden- 
tified in the homology search described above) and random- 
ly selected 170 ESTs (78 ESTs from the Fest_F cDNA 
clones, 92 ESTs from the Fest L cDNA clones) among 863 
cDNA clones were designed by using PRIMER3 software 
(Supplemental Table 1). PCR was performed with candidate 
gene and EST primers and genomic DNA isolated from 
KYU and CAN as the template. The amplified products 
were then sequenced to check whether the target region had 
been amplified. PCR was performed under the following 
conditions: 30 ng of genomic DNA was amplified by PCR in 
a 30 ul final reaction volume containing 6 uM each of the 
forward and reverse primers, 6 mM each of deoxynucleotide 
triphosphate (dNTPs), lOx PCR buffer, and 0.75 units of 
Blend Taq DNA polymerase (TOYOBO Co., Ltd., Tokyo, 
Japan). PCR was performed by using the Mastercycler gra- 
dient PCR system (Eppendorf Scientific Inc., Hamburg, 
Germany). The PCR program consisted of a preliminary de- 
naturation step of 2 min at 94 °C; followed by 40 cycles of 
30 s at 94°C (denaturation), 30 s at gradient temperatures 
(50 to 70°C) (annealing), and 2 min at 72°C (extension); 
with a final extension for 7 min at 72°C. PCR products 
were analyzed by electrophoresis in 8% acrylamide (29 : 1 
acrylamide:bis-acrylamide) gels. For sequencing, the PCR 
products within a single band were purified by using the 
QIAquick PCR purification kit (Qiagen Inc.) and then se- 
quenced by using the BigDye Terminator ver. 3.1 Cycle Se- 
quencing Kit (Life Technologies Corp.) and an ABI PRISM 



3100 Genetic Analyzer (Life Technologies Corp.). 

Primer sets that produced large differences in amplicon 
length (detected on 8% acrylamide gels) between parental 
lines were developed as EST markers. In cases where no dif- 
ference in amplicon length could be detected, we searched 
for single-nucleotide polymorphisms (SNPs) by sequencing, 
and developed cleaved amplified polymorphic sequence 
(CAPS) markers in regions where SNPs were found within 
restriction-enzyme recognition sites. When no SNP could be 
found within the recognition sequences of restriction en- 
zymes, derived cleaved amplified polymorphic sequence 
(dCAPS) markers were developed instead. The SNP searches 
were performed by using Phred/Phrap/Consed software 
(Gordon and Green 1998). The selection of restriction en- 
zymes to detect SNPs and design of mismatch primers for 
dCAPS was conducted with dCAPS Finder 2.0 (Neff et al. 
2002; http://helix.wustl.edu/dcaps/dcaps.html). 

To detect polymorphisms of CAPS and dCAPS markers, 
PCR products produced by using EST primers or dCAPS 
primers were digested with a corresponding restriction en- 
zyme in 12 reaction volumes containing 10 ul of the PCR 
products, 6 units restriction enzyme, and 1 .2 ul of the sup- 
plied lOx buffer. The reaction temperature was set to the 
recommended temperature for each restriction enzyme, and 
the reaction time was 3 h. PCR and electrophoresis condi- 
tions were the same as described above. 

The candidate gene and EST markers were assayed in the 
F 4 mapping population. A framework linkage map for QTL 
analysis was constructed with the software JOINMAP 4.0 
(Van Ooijen 2006) by using genotype data for each marker. 
The fit to the Mendelian segregation ratio for each marker 
was tested with the Chi-square test (P < 0.05). Markers were 
assigned to linkage groups with a logarithm of odds (LOD) 
threshold of 3.0. 

QTL analysis 

The F 4 population was cultivated in a growth cabinet and 
scored for photoperiod sensitivity under a 15.5 h photoperi- 
od, as described above. QTL analysis was performed by us- 
ing the software MAPQTL version 5 (Van Ooijen 2005). In- 
terval mapping analysis was performed (Lander and Botstein 
1989) to locate preliminary QTL positions on the map. Inter- 
val mapping was used to select markers significantly associ- 
ated with the trait to constitute an initial set of cofactors. A 
backward-elimination procedure was applied to the initial 
set of cofactors. Only significant markers (P < 0.02) were 
used as cofactors in the multiple QTL method (MQM) 
(Jansen 1993, Jansen and Stam 1994) analysis for QTL de- 
tection. After the selection of cofactors, MQM analysis was 
performed. A 1000-permutation test was applied to each 
data set to establish the LOD value to be used for deter- 
mining the significance (7 3 <0.05) of identified QTLs 
(Churchill and Doerge 1994). The phenotypic variance ex- 
plained by a single QTL was calculated as the square of the 
partial correlation coefficient (r 2 ) with the observed vari- 
able, adjusted for cofactors. When candidate gene markers 
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could not be mapped, the association between each of candi- 
date gene markers and photoperiod sensitivity was analyzed 
with ANOVA. To evaluate digenic interactions between 
pairs of QTLs, we performed two-way ANOVA for two- 
locus QTL genotypes, as represented by the genotypes of 
markers nearest QTLs These ANOVA were calculated by 
using the program JMP 6.0. 

Results 
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Natural variation in photoperiod sensitivity between autog- 
amous lines 

The frequency distributions of days-to-flowering of KYU 
and CAN are shown in Fig. 1. Under the natural long-day 
condition, the days-to-flowering distributions of the KYU 
(range 22-38 days) and CAN (21-32 days) populations were 
similar (Fig. 1 A), and the means (24.7 and 24.9 for KYU and 
CAN, respectively) were not significantly different 
(P = 0.681) (Fig. 1A). On the other hand, under the artificial 
long-day condition with 1 5.5 h day length, the mean days- 
to-flowering in KYU (mean 24.6; range 24-31 days) was 
significantly shorter (P < 0.0001) than that of CAN (mean 
48.0; range 41-59 days) (Fig. IB). 

Search for photoperiod-sensitivity candidate gene regions 

In tBLASTX searches against the reference sequences of 
proteins in Arabidopsis, three Fest_L EST regions 
(Fest_L0268, Fest_L0327 and Fest_L0352) showed high ho- 
mology with photoperiod-sensitivity genes in Arabidopsis. 
Fest_L0268 showed high homology with CONSTANS-LIKE 
3 (COLS), CONSTANS-LIKE I (COLL) and CONSTANS- 
LIKE 2 (COL2); Fest_L0327, with CIRCADIAN CLOCK 
ASSOCIATED I (CCA1) and LATE ELONGATED HYPO- 
COTYL (LHY); and Fest_L0352, with EARLY FLOWERING 
3 (ELF 3) (Table 1). Among the sequences with homology to 
multiple genes, Fest_L0268 showed the highest homology 
with COL3, and Fest_L0327 showed the highest homology 
with CCAI. We considered these three EST regions as can- 
didate genes related to photoperiod sensitivity in common 
buckwheat, and designated them FeCOL3, FeCCAl, and 
FeELF3, respectively. None of the Fest_F regions showed 
high homology to photoperiod genes in Arabidopsis. 

Development of linkage map using EST markers 

We developed primers for the three candidate genes 
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Fig. 1. Variation in days-to-flowering between the parental lines. (A) 
Variation in days-to-flowering of parental lines under natural long-day 
condition (-14.5 h day length). (B) Variation in days-to-flowering of 
parental lines under artificial long-day condition (15.5hday length). 
KYU, 02AL1 13(Kyukei SC2)LH.self, CAN, C0408-0 RP; P value in 
each figure showed significance of difference by t-test. 

(FeCOL3, FeCCAl and FeELFS) and randomly selected 
170 ESTs (78 from the Fest_F group, 92 from the Fest_L 
group) from among the original 863 EST regions (Supple- 
mental Tables 1 and 2 showed results of BLASTX for ran- 
domly selected 170 EST regions). Among these, we con- 
firmed polymorphism between KYU and CAN for the three 
candidate gene and 115 ESTs regions. Sixteen of the 115 
polymorphic EST regions included large indels, whereas the 
three candidate gene and 99 ESTs regions included small in- 
dels and SNPs. The 16 EST regions containing large indels 
were used as length markers; in addition, we developed three 
candidate gene and 60 SNP-based markers (39 CAPS and 21 
dCAPS markers) from among the 99 ESTs markers. Infor- 
mation on the three candidate gene markers and 76 EST 
markers is shown in Table 2. 

Through the linkage mapping of these markers, the posi- 
tions of 63 EST markers were estimated. Nine other EST 
markers showed departure from Mendelian segregation ra- 
tios, and four EST markers could not be mapped (Table 2 
and Fig. 2). The map consisted of 9 linkage groups (LG) 



Table 1. Confirmed sequences showing high homology with photoperiod-sensitivity genes in Arabidopsis 



EST region 


Gene product 


Gene ID 


Bit score 


E-value 


Fest_ 


L0268 


COL3 (CONSTANS-LIKE 3); protein binding / transcription factor/ zinc ion binding 


817016 COL3 


124 


4.00E-29 






COL1 (constans-like I); transcription factor/ zinc ion binding 


831442 COL1 


112 


2.00E-25 






COL2 (constans-like 2); transcription factor/ zinc ion binding 


821298 COL2 


112 


2.00E-25 


Fest_ 


L0327 


CCAI (CIRCADIAN CLOCK ASSOCIATED 1); DNA binding / transcription activator/ 


819296 CCAI 


142 


4.00E-34 






transcription factor/ transcription repressor 












LHY (I ATE ELONGATED HYPOCOTYL); DNA binding / transcription factor 


839341 LHY 


137 


5.00E-33 


Fest_ 


L0352 


EIF3 (EARLY FLOWERING 5); protein C-terminus binding / transcription factor 


817134 ELF3 


104 


1.00E-22 
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Fig. 2. Linkage map of common buckwheat. 
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ranging from 2.6 cM (LG2) to 87.6 cM (LG3) and covered 
31 1.6 cM in total. The segregation ratios for the three candi- 
date gene markers were not significantly different from the 
expected Mendelian ratios. FeCOL3 and FeCCAl markers 
were mapped to LG9 and LG6, respectively, but FeELF3 
could not be mapped (Table 2 and Fig. 2). 

QTL analysis for photoperiod sensitivity 

The segregation of photoperiod sensitivity in the F 4 prog- 
enies was observed as a wide but unimodal distribution. 
Transgressive segregation in the direction of late flowering 
was observed (Fig. 3). 

When we performed interval mapping, two regions were 
associated with photoperiod sensitivity: one in LG3 (25.3 to 
50.3 cM) and one in LG7 (0 to 15.1 cM) (LOD score >2.4, 
P< 0.05). The results of cofactor analysis indicated a signif- 
icant difference (P < 0.02) at eight markers (Fest_L0490_2 
and Fest_L0472_l in LG1, Fest_L0606_4 in LG3, 
Fest_F0077_2 and Fest LOlOl l in LG4, Fest_L0337_6 in 
LG7 and Fest_L0242_l and Fest_L0053_3 in LG8). In the 
MQM analysis, two QTLs (Fest_L0606_4 and 
Fest_L0337_6 in LG3 and 7, respectively) showed relatively 
high LOD values (LOD score >2.4, P < 0.05) (Table 3 and 
Fig. 4). For each of these two QTLs, the KYU alleles had 
negative additive effects on days-to-flowering (-8.2 days 



30 

» 25 

5-15 
o 

6 10 



KYU 
24.6±2.02 



CAN 
48.0±5.50 



n,n.n., 



16- 21- 26- 31- 36- 41- 



46- 51- 56- 61- 66- 71- 
Days-to-flowering 



76- 81- 86- 91- 96- 



Fig. 3. Distribution of days-to-flowering in F 4 progeny (n = 1 00) of a 
cross of 02AL113(Kvukei SC2)LH.self (KYU) x C0408-0 RP (CAN). 
Arrows indicate means ± SD for each parent, n.f, non-flowering. 



and -6.8 days, for Fest_L0606_4 and Fest_L0337_6, re- 
spectively) and explained about 20.0% and 14.2% of the 
phenotypic variation in the F 4 progenies, respectively. 

In the candidate gene markers, FeCCAl and FeCOL3 
showed low LOD values (LOD score 0.13 and 0.1 1, respec- 
tively) as a result of the MQM analysis (Table 3). FeELF3 
could not be mapped. Thus, we searched for associations be- 
tween FeELF3 and photoperiod sensitivity by using the 
ANOVA. Table 4 shows the result of the ANOVA for asso- 
ciation between the FeELF3 and photoperiod sensitivity. It 
was shown that FeELF3 was significantly related to photo- 
period sensitivity (P- 0.028). For FeELF3 markers, CAN 
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Table 3. Quantitative trait loci (QTLs) detected for photoperiod sensi- 
tivity in the 02AL113(Kyukei SC2)LH.self (KYU) x C 040 8-0 RP 
(CAN) F4 population by multiple QTL analysis 



LG 


Position 
(cM) 


Locus 




LOD 


Additive 
effect 


PVE (%) 


3 


34.0 


Fest_L0606 


4 


5.2 


-8.2 


20.0 


7 


0.0 


Fest_L0337_ 


6 


3.8 


-6.8 


14.2 


6 


21.5 


FeCCAl 




0.13 


1.2 


0.4 


9 


0.6 


FeCOL3 




0.11 


-1.1 


0.4 


U 




FeELF3 











U, unknown genome location 

Table 4. Quantitative trait loci (QTLs) detected for photoperiod sensi- 
tivity in the 02AL113(Kyukei SC2)LH.self (KYU) x C 040 8-0 RP 
(CAN) F 4 population by ANOVA. KYU and CAN indicate average of 
days-to-flowering for the 02AL1 13(Kyukei SC2)LH.self allele and 
homozygosity for the C0408-0 RP allele, respectively. 



LG 


Position 
(cM) 


Locus 




F value 


P value 


Average of days- 
to-flowering 












KYU 


CAN 


U 




FeELF3 




3.72 


0.028 


51.5 


42.9 


3 


34.0 


Fest_L0606 


4 


10.40 


<0.001 


39.3 


54.6 


7 


0.0 


Fest_L0337 


6 


6.17 


0.003 


40.6 


53.1 



U, unknown genome location 



allele showed smaller average of days-to-flowering than 
KYU allele, the difference was 8.6 days while KYU allele 
had early average of days-to-flowering (the difference was 
15.3 and 12.5 days, for Fest_L0606_4 and Fest_L0337_6, 
respectively) in the two QTL's. 

To test for the existence of QTL interactions, genotype 
data for Fest_L0606_4, Fest_L0337_6 and FeELF3 were 
used in a two-way ANOVA to test the phenotypic difference 
among the four homozygous genotype classes (Fig. 5). A 
significant interaction was detected between Fest_L0606_4 
and Fest_L0337_6 LP < 0.01). On the other hand, interac- 
tions between Fest_L0606_4 and FeELF3 and between 
Fest_L0337_6 and FeELF3, were not significant (,P = 0.81 
and P = 0.89, respectively). 
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Fig. 4. LOD score for QTLs detected in multiple QTL analysis. The 
QTL likelihood maps for each linkage group were obtained by using 
the MQM procedure of MapQTL, fixing one marker cofactor per puta- 
tive QTL. Selected markers used as cofactors in analysis are indicated 
by black rectangles along the x-axes. The linkage group number is in- 
dicated at the top of each graph. Horizontal dashed lines indicate the 
2.4 (P < 0.05) LOD score threshold. 



Discussion 

A number of studies have been conducted for investigating 
photoperiod sensitivity in common buckwheat. Minami 
(1985) suggested that photoperiod sensitivity of common 
buckwheat is controlled by the action of multiple genes. 
However, neither genes nor candidate regions related to 
photoperiod sensitivity had been identified until now. In the 
present study, a difference in days-to-flowering between the 
KYU and CAN populations was expressed under a 15.5 h 
photoperiod, but no difference was observed under the 14.5 h 
photoperiod. These results showed that the difference in 
days-to-flowering between KYU and CAN observed under 
the 15.5 h photoperiod was caused by a difference of photo- 



period sensitivity: the KYU population had weaker photo- 
period sensitivity than the CAN population. The KYU line was 
bred by interspecific hybridization among F. homotropicum 
(an autogamous species) and F. esculentum cv. Botansoba 
which is cultivated in high-latitude regions (Matsui 2006); 
similarly, the CAN line was also bred by crossing 
F. homotropicum and F. esculentum which is cultivated in 
high-latitude regions (Wang and Campbell 1 998). Thus, we 
assumed that both lines had been selected under long-day 
conditions. This may be why the difference in photoperiod 
sensitivity between KYU and CAN was not evident unless 
an artificially long day length (15.5 h) was applied. 

In Arabidopsis, the photoperiod-sensitivity pathway 
consists of photoreceptors, the circadian clocks, and floral 
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Fig. 5. Mean values for days-to-flowering in different combinations of genotype classes for Fest_L0606_4 and Fest_L0337_6. The x-axis indi- 
cates the genotype at Fest_L0606-4 and Fest_L0337_6: KYU and CAN indicate homozygosity for the 02AL1 13(Kyukei SC2)LH.self allele and 
homozygosity for the C0408-0 RP allele, respectively. The genotype of Fest_L0337_6 and FeELF3 are indicated by circles (homozygous for 
KYU allele) or squares (homozygous for CAN allele). P value in each figure showed results of significance for interaction by two way ANOVA. 



promoters (Putterill et al. 2004). In many species with photo- 
period sensitivity, such as rice, wheat (Triticum aestivum), 
and morning glory (Pharbitis nil), genes orthologous to 
photoperiod-sensitivity genes found in Arabidopsis were 
confirmed to be related to heading date or flowering time 
(Hayama et al. 2002, 2003, Kojima et al. 2002, Liu et al. 
2001a, Nemoto et al. 2003, Yano et al. 2000). In the present 
study, we found three such candidate gene regions 
(FeCOL3, FeCCAl and FeELF3) among 863 cDNA clones. 
In Arabidopsis, CCA1 and ELF3 are genes that connect the 
circadian clock to the photoperiod-sensitivity pathway 
(Covington et al. 2001, Hicks et al. 2001, Liu et al. 2001b, 
Wang et al. 1997, Wang and Tobin 1998), and COL3 is from 
a CO-related gene family that connects floral promoters to the 
photoperiod-sensitivity pathway (Putterill et al. 1995). Here, 
we found that the FeELF3 candidate gene region in common 
buckwheat was strongly associated with photoperiod sensi- 
tivity. Although the candidate genes identified in the present 
study would only be part of the total group of genes included 
in a photoperiod-sensitivity pathway, our report is the first in 
which candidate genes identified based on their role in photo- 
period sensitivity in another species were associated with 
photoperiod sensitivity in common buckwheat. The results 
indicate that it would be efficient to study the photoperiod- 
sensitivity pathway by using other candidate genes and ana- 
lyzing their effects on photoperiod sensitivity. Although mu- 
tations in intron region were reported to be associated with 
gene expression (Fiume et al. 2004, Fu et al. 2005, Isshiki 
et al. 1998, Samadder et al. 2008), it was not clear whether 
FeELF3 marker based on SNP in an intron region directly 
associated with photoperiod sensitivity right now. If SNP's 
associated with photoperiod sensitivity in direct will be de- 
tected in the future, the allelic variations in FeELF3 gene 
region will become more informative information. 

In the MQM analysis, we detected two QTLs, 
Fest_L0606_4 and Fest_L0337_6, which explained about 
20.0% and 14.2% the variation observed in the F 4 progenies, 
respectively. In the present study, we constructed a linkage 
map for QTL analysis consisting of 9 linkage groups con- 



taining 63 ESTs and two candidate gene regions and cover- 
ing 311.6 cM in total. This map was the first linkage map to 
include ESTs and candidate gene regions thought to be relat- 
ed to photoperiod sensitivity in common buckwheat. Unfor- 
tunately, we assumed that the genome coverage of this map 
is not so high because Konishi and Ohnishi (2006) con- 
structed a linkage map of buckwheat covering about 900 cM 
by using AFLP and SSR markers. Thus, result of MQM 
analysis at unmapped FeELF3 candidate gene region was not 
shown. If the genome coverage of this map will be higher, 
other QTLs included FeELF3 candidate gene region might 
be detected. 

The genetic analysis of quantitative traits with molecular 
markers makes it possible to evaluate epistatic interactions 
between individual loci. A digenic interaction between the 
QTLs detected in the MQM analysis was observed in this 
study. The additive effect of the KYU allele of 
Fest_L0337_6 was estimated to be -6.8 days. However, the 
observed two-locus interaction between Fest_L0606_4 and 
Fest_L0337_6 suggests that the effects of Fest_L0337_6 are 
changed depending on the allele present at Fest_L0606_4. In 
Arabidopsis and rice, interactions between genes related to 
photoperiod sensitivity have been reported (Putterill et al. 
2004, Yano et al. 1997). These results suggested that photo- 
period sensitivity in common buckwheat was controlled by 
both main and interaction effects of the QTLs detected. 
Thus, it is necessary to analyze both the main effects and the 
interactions between QTLs to fully understand photoperiod 
sensitivity. 

In summary, our study is the first report that photoperiod 
sensitivity in common buckwheat is controlled by the action 
of multiple genes, including interaction effects between 
these genes. The photoperiod sensitivity was controlled by at 
least three loci, including candidate genes identified as being 
part of the photoperiod-sensitivity pathway in Arabidopsis. 
These findings will serve as a stepping-stone to future stud- 
ies on the role of photoperiod sensitivity in local adaptation 
and ecological breeding in common buckwheat. 
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